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Abstract 


The  CTH  code  was  used  to  simulate  sympathetic  detonation  experiments  conducted  using 
105-mm  projectiles  separated  by  buffers  of  Plexiglas,  rolled  homogeneous  armor  (RHA),  and  mild 
steel  Buffer  thicknesses  near  the  ejq)erimental  sympathetic  detonation  threshold  were  simulated. 
No  propagation  criterion  associated  with  pressure  loading  was  revealed  by  the  results  of  this 
study.  CTH  con:5)utations  were  also  made  to  confirm  the  presence  of  finite  rise-time  or  “ramp” 
waves  in  the  acceptor  explosive  (previously  observed  in  Lagrangian  computations)  with  layered 
buffers  of  Plexiglas  and  steel.  The  CTH  results  show  more  structure  in  the  ramp  wave  and 
shorter  rise  times  than  those  previously  obtained.  The  History  Variable  Reactive  Bum  (HVRB) 
explosive  initiation  model  in  the  CTH  code  was  exercised  in  simulations  of  similar  configurations. 
These  simulations  showed  that  initiation  may  occur  during  propagation  of  the  incident  shock 
wave  through  the  acceptor  or  after  its  reflection  from  the  casing  at  the  back  of  the  acceptor. 
Predicted  buffer  thicknesses  required  to  prevent  sympathetic  detonation  are  much  less  than  those 
determined  in  the  experiments,  indicating  that  mechanisms  in  addition  to  shock  initiation 
contribute  to  sympathetic  detonation. 


TABLE  OF  CONTENTS 


Page 

LISTOFHGURES .  v 

LIST  OF  TABLES .  ^ 

1.  INTRODUCTION .  1 

2.  DESCRIPTION  OF  CTH  .  1 

3.  BRIEF  DESCRIPTION  OF  THE  EXPERIMENTS .  2 

4.  era  SIMULATION  OF  THE  EXPERIMENTS .  2 

5.  COMPARISON  OF  era  AND  STEALTH  COMPUTATIONS  .  4 

6.  HVRB  PREDICTIONS  OF  SYMPATHETIC  DETONATION .  5 

7.  SUMMARY  AND  CONCLUSIONS . 5 

8.  REFERENCES  .  27 

DISTRIBUTION  LIST .  29 

REPORT  DOCUMENTATION  PAGE .  33 


m 


Intentionally  left  blank. 


IV 


LIST  OF  HGURES 


Figure  2age 

1.  Experimental  arrangement  used  by  Boyle .  7 

2.  2-D  CTH  simulation .  9 

3.  76-nmi  Plexiglas  buffer  contour  plot  sequence .  10 

4.  76-mm  Plexiglas  buffer  pressure  history  .  11 

5.  76-mm  Plexiglas  buffer  impulse  history .  12 

6.  Plexiglas  and  mild  steel:  comparison  of  pressure  contours  .  13 

7.  Plexiglas  and  mild  steel:  comparison  of  pressure  and  impulse  histories .  14 

8.  Mild  steel  and  RHA:  comparison  of  pressure  contours  at  the  experimental 

threshold .  1^ 

9.  Mild  steel  and  RHA:  comparison  of  pressure  and  impulse  histories  at  the 

experimental  threshold .  16 

10.  Mild  steel  and  RHA:  comparison  of  pressure  contours  at  the  same  buffer 

thickness  .  17 

1 1 .  Mild  steel  and  RHA:  comparison  of  pressure  and  impulse  histories  at  the  same 

buffer  thickness  .  1^ 

12.  Results  of  the  STEALTH  sympathetic  detonation  simulation  with  layered  buffers  .  20 

13.  Resultsof  the  CTH  sympathetic  detonation  simulation  with  layered  buffers .  21 

14.  Unbuffered  configuration  with  12.7-mm  separation:  pressure  and  reaction 

variable  contours .  22 

15.  6.4-mm  Steel  buffer  in  direct  contact  with  donor  and  acceptor  munitions:  pressure 

and  reaction  variable  contours  .  23 

16.  6.4-mm  Plexiglas  buffer  centered  in  12.7-mm  space:  pressure  and  reaction  variable 

contours .  24 


V 


Intentionally  left  blank. 


VI 


LIST  OF  TABLES 


Table  Cage 

1.  Experimental  Results .  8 

2.  Confuted  Pressure  and  Impulse  at  the  Experimental  Thresholds .  19 

3.  Sympathetic  Detonation  Predictions  for  the  Contact  Configuration .  25 

4.  HVRB  Results  for  12.7-mm  Munition-to-Munition  Separation  Computations  ...  26 


vii 


Intentionally  left  blank. 


1.  INTRODUCTION 


Sympathetic  detonation  of  munitions  occurs  when  a  detonating  munition  (die  donor)  initiates 
detonation  in  a  neighboring  munition  (the  acceptor).  Shock  due  to  casing  impact  is  the  primary 
propagation  mechanism  for  unshielded  munitions  at  close  range,  while  at  longer  range,  fragment 
penetration  is  the  primary  mechanism.  Shielding  (buffers)  between  munitions  reduces  the 
incidence  of  sympathetic  detonation  by  reducing  shock  stimulus  and  impeding  fragments. 

The  present  study  employed  the  CTH  code  (Hertel  et  al.  1993)  and  focused  on  three 
objectives.  The  first  was  to  determine  if  any  similarities  in  pressure-loading  exist  when  using  CTH 
to  simulfltft  buffered  sympathetic  detonation  experiments  with  buffer  thicknesses  near  the 
sympathetic  detonation  threshold,  as  determined  in  experiments.  The  second  objective  was  to 
confirm  the  presence  of  finite  rise-time  or  “ramp”  waves  in  the  acceptor  explosive  previously 
observed  using  the  nonreactive  Lagrangian  code  STEALTH  (Electric  Power  Research  Institute 
1981).  The  third  objective  was  to  exercise  a  CTH  reactive  model  within  the  acceptor  explosive. 

2.  DESCRIPTION  OF  CTH 

The  continuum  mechanics  solver  CTH  is  an  ongoing  project  of  the  Sandia  National 
Laboratory.  It  is  intended  to  provide  capabilities  for  modeling  dynamics  of  multidimensional 
systems  with  multiple  materials,  large  deformations  and  strong  shock  waves.  Finite-difference 
analogs  of  the  Lagrangian  equations  of  momentum  and  energy  conservation  are  employed  with 
continuous  rezoning  to  construct  Eulerian  differencing.  Shock  and  detonation  waves  are  treated 
using  the  method  of  artificial  viscosity.  CTH  makes  use  of  analytic  (Mie-Gruneisen,  JWL,  etc.) 
and  tabular  (Sesame)  equations  of  state  as  well  as  modem  constitutive  models  (Johnson-Cook, 
Zerilli-Armstrong)  including  fracture  (void  insertion).  In  addition,  there  are  three  reaction  models 
(Programmed  Bum,  CJ  Volume  Bum  and  History  Variable  Reactive  Bum  [HVRB])  and  two 
porosity  models.  A  variety  of  plots  of  computed  results  can  be  produced.  These  options  provide 
an  opportunity  to  treat  complex  material  behavior  including  melting,  vaporization,  solid  phase 
transitions,  chemical  reaction,  and  electronic  excitation  and  ionization. 
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3.  BRIEF  DESCRIPTION  OF  THE  EXPERIMENTS 


The  experimental  arrangement  (Boyle  1995)  is  illustrated  in  Figure  1.  Both  the  donor  and  the 
acceptor  munitions  were  105-mm  projectiles.  In  each  experiment,  the  donor  was  filled  with 
Composition  B  and  the  acceptor  with  pentolite.  Each  donor  was  initiated  at  a  point  on  the  side 
farthest  from  the  acceptor,  halfway  between  its  top  and  bottom.  Potted  carbon  gauges  and  steel 
witness  plates  were  used  to  record  the  results.  Buffer  materials  included  polyethylene,  mild  steel, 
and  rolled  homogeneous  armor  (RHA).  All  buffers  were  203  mm  wide  and  510  mm  high.  The 
thickness  of  the  buffers  was  varied  in  order  to  determine  a  sympathetic  detonation  threshold. 

The  experimental  results  are  summarized  in  Table  1.  The  “GOs”  represent  configurations  in 
which  sympathetic  detonation  of  the  acceptor  occurred,  while  the  “NOGOs”  indicate  that 
sympathetic  detonation  did  not  occur.  Threshold  ranges  obtained  from  these  results  determined 
the  buffer  thicknesses  we  simulated  computationally. 

4.  era  SIMULATION  OF  THE  EXPERIMENTS 

The  2-D  era  simulation  (representing  the  cross-sectional  plane  in  which  the  donor  is 
initiated)  is  illustrated  in  Figure  2.  The  primary  effect  of  simulating  a  3-D  experiment  in  2-D  is  to 
decrease  the  rate  of  decay  of  spherically  divergent  waves  to  that  of  cylindrically  divergent  waves. 
This  produces  somewhat  higher  stimulus  levels  in  the  2-D  case.  The  experiments  were  simulated 
with  buffer  thicknesses  near  the  sympathetic  detonation  threshold  for  each  buffer  material.  The 
dimensions  were  consistent  with  Ml  105-mm  projectiles.  Casings  were  lO-mm-thick  mild  steel. 
The  donor  munition  was  filled  with  Composition  B,  modeled  using  a  JWL  equation  of  state  for 
detonation  products.  The  acceptor  fill  was  "inert"  pentolite,  represented  by  a  JWL  equation  of 
state  with  estimated  constants.  The  donor  munition  was  initiated  at  the  side  farthest  firom  the 
acceptor  using  the  programmed  bum  model.  All  buffers  were  203  mm  wide.  Plexiglas  was 
substituted  for  polyethylene,  which  was  used  as  a  buffer  material  in  some  of  the  original 
experiments.  The  buffer  thicknesses  were  70  mm  and  76  mm  for  the  Plexiglas  buffers,  32  mm  and 
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38  mm  for  the  mild-steel  buffers,  and  38  mm  and  5 1  mm  for  the  RHA  buffers.  In  addition  to  the 
threshold  computations,  a  computation  was  made  with  a  38-mm-thick  RHA  buffer  to  enable 
direct  comparison  with  the  mild-steel  computation  at  the  same  thickness. 

Figure  3  shows  a  sequence  of  pressure  contour  plots  for  the  Plexiglas  buffer.  The 
programmed  bum  detonation  can  be  seen  propagating  across  the  donor  munition  in  the  first 
frame,  while  subsequent  frames  Ulustrate  a  shock  wave  propagating  through  the  buffer  and  casing 
of  the  acceptor  and  into  the  acceptor  explosive.  The  buffer  is  nearly  obliterated  and  the  acceptor 
is  significantly  deformed  at  late  times.  A  crack  can  be  seen  in  the  casing  at  the  back  of  the 
acceptor  in  the  last  firame.  Figures  4  and  5  illustrate  pressure  and  impulse  history  plots, 
respectively,  from  the  same  computation,  taken  at  a  point  just  inside  die  acceptor  casing.  The 
loading  consists  of  an  initial  pulse  about  40  ps  wide,  followed  by  significant  late-time 
reverberations.  The  impulse  history  plot  indicates  that  about  half  the  ultimate  impulse  delivered  is 
associated  with  the  reverberations. 

Figure  6  is  a  comparison  of  pressure  contour  plots  for  the  Plexiglas  buffer  with  sumlar  plots 
for  a  mild-steel  buffer.  The  shock  arrives  at  the  acceptor  explosive  15  ps  sooner  in  die  mild-steel 
buffer  simulation,  due  to  the  smaller  thickness  of  and  higher  wave  speed  in  the  steel.  Otherwise, 
the  results  from  both  the  Plexiglas  and  mild-steel  buffer  simulations  appear  similar.  At  late  times, 
some  cracking  on  the  front  of  the  acceptor  near  the  center  can  be  observed  with  the  steel  buffer 
and  acceptor  casing,  but  no  cracking  is  seen  at  the  back  of  the  casing.  As  illustrated  in  Figure  7, 
significant  differences  in  pressure-loading  are  indicated  in  the  pressure  and  impulse  history  plots. 
With  the  mild-steel  buffer,  the  initial  pulse  is  higher  and  narrower  and  the  late-time  pressures 
remain  insignificant  Thus,  though  the  impulse  associated  with  the  initial  pulse  is  higher  with  the 
mild-steel  buffer,  the  ultimate  impulse  is  higher  with  the  Plexiglas  buffer. 

Figure  8  illustrates  a  comparison  of  pressure  contour  plots  for  the  38-mm  mild-steel  buffer 
and  the  5 1  -mm  RHA  buffer.  With  the  RHA  buffer,  a  more  definite  cracking  response  can  be  seen 
in  the  acceptor  casing  on  the  side  near  the  buffer,  but  otherwise  the  results  appear  similar.  A 
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comparison  of  the  pressure  and  impulse  history  plots  of  the  RHA  and  mild-steel  buffers,  as  shown 
in  Figure  9,  reveals  no  significant  differences  except  that  the  pressures  are  somewhat  higher  and 
die  ultimate  impulse  is  significandy  higher  with  the  mild-steel  buffer. 

Figure  10  shows  a  comparison  of  pressure  contour  plots  between  a  38-mm  RHA  buffer  and  a 
38-mm  mild-steel  buffer.  With  the  buffer  thicknesses  equal,  the  RHA  shows  greater  resistance  to 
deformation  and  cracks  clearly  at  late  times.  The  pressure  and  impulse  history  plots  are  similar, 
as  shown  in  Figure  11.  The  initial  peak  pressure  is  slighdy  higher  and  the  ultimate  impulse  is 
somewhat  lower  with  the  RHA  buffer. 

The  pressure  and  impulse  results  at  the  experimental  threshold  are  summarized  in  Table  2. 

The  impulse  at  the  end  of  the  first  pulse,  as  well  as  the  ultimate  impulse  (defined  as  the  value  at 
250  ps)  is  included.  As  shown  in  the  table,  the  peak  pressure  is  the  same  for  Plexiglas  and  RHA, 
but  higher  for  mild  steel.  The  first-pulse  impulse  is  different  for  all  three  buffer  materials.  The 
ultimate  impulse  values  are  similar  for  Plexiglas  and  mild  steel,  but  lower  for  RHA. 

5.  COMPARISON  OF  CTH  AND  STEALTH  COMPUTATIONS 

CTH  computations  were  performed  to  compare  with  the  results  obtained  in  a  previous 
sympathetic  detonation  study,  which  used  the  nonreactive  Lagrangian  code  STEALTH 
(Starkenberg  et  al.  1987).  In  the  STEALTH  study,  ramp  waves  as  long  as  25  ps  were  observed 
in  the  acceptor  explosive  when  layered  buffers  were  used  between  munitions.  In  one  case,  the 
layered  buffers  consisted  of  10  mm  of  Plexiglas  sandwiched  between  two  10-nun  layers  of  steel. 
Figures  12  and  13  show  pressure  waveforms  produced  in  the  acceptor  munition  in  STEALTH  and 
CTH  simulations  of  the  same  configmation.  The  results  obtained  witii  CTH  show  more  stmcture 
in  the  ramp  wave  and  shorter  rise-times  than  those  obtained  with  STEALTH. 
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6.  HVRB  PREDICTIONS  OF  SYMPATHETIC  DETONATION 


The  HVRB  explosive  initiation  model  (Kerley  1992)  was  exercised  in  two  series  of  CTH 
simulations.  In  one  series,  the  donor  and  acceptor  were  in  direct  contact  with  the  buffer.  In  the 
other,  the  donor  and  acceptor  were  separated  by  a  fixed  12.7-mm  space,  and  the  buffer  was 
centered  in  that  space.  In  both  series,  the  acceptors  were  filled  with  reactive  Composition  B,  and 
buffo"  thicknesses  were  decreased  imtil  sympathetic  detonation  occurred. 

% 

Figure  14  shows  a  sequence  of  pressure  and  reaction  variable  contour  plots  for  imbuffered 
projectiles  separated  by  a  12.7-mm  space.  Pressure  contours  are  shown  on  the  right  side,  and 
reaction  variable  contours  are  shown  on  the  left  side  of  each  plot.  In  this  case,  initiation  of  the 
acceptor  occurs  as  the  incident  shock  wave  propagates  through  it. 

A  sequence  of  contour  plots  for  a  6.4-mm-thick  steel  buffer  in  direct  contact  with  the  donor 
and  acceptor  munitions  is  shown  in  Figure  15.  This  illustrates  a  configuration  that  does  not 
produce  sympathetic  detonation. 

Figure  16  shows  a  sequence  of  pressure  and  reaction  variable  contour  plots  for  a  6.4-mm 
Plexiglas  buffer  centered  in  a  12.7-mm  space.  Initiation  of  the  acceptor  occurs  after  reflection  of 
the  shock  wave  at  55  ps. 

The  results  obtained  in  the  HVRB  simulations  are  summarized  in  Tables  3  and  4.  In  all  cases, 
the  buffer  thicknesses  required  to  prevent  sympathetic  detonation  are  less  than  those  required  in 
the  experiments. 

7.  SUMMARY  AND  CONCLUSIONS 

We  simulated  sypathetic  detonation  experiments  with  buffers  at  the  threshold  thickness.  No 
propagation  critericm  associated  with  pressure-loading  was  revealed  by  the  results  of  this  study. 
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Buffers  may  prevent  sympathetic  detonation  by  influencing  other  aspects  of  the  loading  in  addition 
to  pressure.  We  observed  that  low-impedance  Plexiglas  buffers  tend  to  spread  out  the  pressure¬ 
loading  over  a  longer  period  of  time,  while  the  high-strength  steel  and  RHA  buffers  reduce 
deformation  and  damage. 

When  layered  buffers  are  used,  the  results  obtained  with  CTH  show  more  structure  in  the 
ramp  wave  and  shorter  rise-times  than  those  previously  obtained  with  STEALTH.  The  response 
of  explosives  to  ramp  waves  is  not  well  understood,  but  it  has  been  observed  to  be  less  violent 
than  to  a  shock  wave  of  the  same  strength. 

Simulations  with  the  HVRB  explosive  initiation  model  show  that  shock  initiation  may  occur 
during  propagation  of  the  incident  shock  wave  through  the  acceptor  or  after  its  reflection  from 
the  casing  at  the  back  of  the  acceptor.  Buffers  required  to  prevent  sympathetic  detonation  in  the 
computations  were  much  thinner  than  those  required  in  the  experiments.  This  is  another 
indication  that  mechanisms  in  addition  to  shock  initiation  contribute  to  S5mipathetic  detonation. 
These  results  contrast  sharply  with  those  obtained  in  an  earlier  study  using  the  Forest  Fire 
initiation  model  and  the  2DE  code  (Starkenberg,  Huang,  and  Arbuckle  1984).  In  that  study,  no 
buffer  that  could  prevent  sympathetic  detonation  of  Composite  B-filled  105-nim  projectiles  was 
found. 
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Table  1.  Experimental  Results 


MATERIAL 

IHICKNESS 

(mm) 

RESULT 

MILD  STEEL 

<38 

4GOs 

38 

NOGO 

38 

NOGO 

>38 

3NOGOS 

RHA 

<51 

2GOs 

51 

GO 

51 

NOGO 

>51 

3NOGOS 

POLYETHYLENE 

<70 

5GOs 

* 

70 

GO 

70 

NOGO 

76 

GO 

76 

NOGO 

76 

NOGO 

* 

>76 

3  NOGOs 

8 


(mm) 


(mm) 


Figure  2.  2-D  CTH  Simulation. 
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PRESSURE  (10®D/CM2) 


Figured  76-mm  Plexiglas  buffer  pressure  history. 
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IMPULSE  {10^G/CM-S) 


HISTORY  POINT  2 


Figures.  76-mm  Plexiglas  buffer  impulse  history. 
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38-mm  Mild  Steel 
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Figure?.  Plexiglas  and  mild  steel:  comparison  of  pressure  and  impulse  histories. 
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MPULSC  (10*C/CM-S)  PRESSURC  (10*D/CM») 


38-mm  Mild  Steel 


51-mm  RHA 


Figure  9.  Mild  steel  and  RHA:  comparison  of  pressure  and  impulse  histories  at  the  experimental 
threshold. 


16 


4  0 

35 

30 

25 

f  20 

V 

15 

10 

5 

0 

40 

35 

30 

25 

1  20 

>“ 

15 

10 

5 

0 

40 

35 

30 

25 

J.  20 
>- 

15 

10 

5 

0 


Figure  10 


impw.se  (io*g/cm-s)  pressure  (t0*0/CM>) 


8 


Table  2.  Computed  Pressure  and  Impulse  at  the  Experimental  Thresholds 


BIJEFER 

BUFFER 

PEAK 

FIRST-PULSE 

ULTIMATE 

MATERIAL 

TfflCKNESS 

PRESSURE 

IMPULSE 

IMPULSE 

(mm) 

(GPa) 

(kPa»s) 

(kPa»s) 

PLEXIGLAS 

70 

0.94 

17.6 

36.0 

76 

16.9 

34.2 

MILD 

32 

1.22 

24.7 

35.1 

STEEL 

38 

1.28 

20.3 

28.8 

RHA 

51 

0.94 

12.9 

16.8 

19 


O.oo  5.00 


0.0 


5.0  20.0  25.0 


Figure  12.  Results  of  the  STEALTH  sympathetic  detonation  simulation  with  layered  buffers. 
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Figure  14.  Unbuffered  configuration  with  12.7-mm  separation:  pressure  and  reaction  variable 
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Table  3.  Sympathetic  Detonation  Predictions  for  the  Contact  Configurations 


BUFFER 

THICKNESS 

(mm) 

STEEL 

PLEXIGLAS 

AIR 

3.2 

GO 

(reflected 
wave,  ~45ps) 

GO 

(reflected 
wave,  ~50ps) 

’  6.4 

NOGO 

NOGO 

GO 

(incident 
wave,  ~23ps) 

12.7 

NOGO 

GO 

(incident 
wave,  ~27ps) 

25 


Table  4.  HVRB  Results  for  12.7 -mm  Munition-to-munition  Separation  Computations. 


BUb'FHR 

TfflCKNESS 

(mm) 

STEEL 

PLEXIGLAS 

AIR 

1.6 

GO 

(incident 
wave,  ~30ps) 

GO 

(incident 
wave,  ~25ps) 

3.2 

NOGO 

GO 

(incident 
wave,  ~25ps) 

6.4 

NOGO 

GO 

(reflected 
wave,  ~55ps) 

12.7 

GO 

(incident 
wave,  ~27us) 
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